Effect of edge decoration on the energy spectrum of semi-infinite lattices 
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Analytical studies ol the effect of edge decoration on the energy spectrum of semi-infinite one- 
dimensional (ID) lattice chain with Peierls phase transition and zigzag edged graphene (ZEG) are 
presented by means of transfer matrix method, in the frame of which the sufficient and necessary 
conditions for the existence of the edge states are determined. For ID lattice chain, the zero-energy 
edge state exists when Peierls phase transition happens regardless whether the decoration exists or 
not, while the non-zero-energy edge states can be induced and manipulated through adjusting the 
edge decoration. On the other hand, the semi-infinite ZEG model with nearest- neighbor interaction 
can be mapped into the ID lattice chain case. The non-zero-energy edge states can be induced 
by the decoration as well, and we can obtain the condition of the decoration on the edge for the 
existence of the novel edge states. 
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I. INTRODUCTION 

One of the interesting phenomena in condensed mat- 
ter physics is the existence of edge states, the proper- 
ties of which are distinct from those of the bulk states. 
There are some examples showing that the system is in- 
sulated in the bulk, while conduction can be allowed by 
edge states on the boundary. The most prominent ones 
are the quantum Hall effect (QHE) 1-4 and the quantum 
spin Hall effect (QSHE) 5 ~ 7 , where the quantization of 
a Hall conductance is tightly associated with the edge 
states^—. 

With the development of graphen o 14 ' 15 in recent years, 
the edge states of graphene models attract much at- 
tention. From the topological view, edge states can 
be induced in the system with different structures 
of edge a 19 i 20 ; in experiments, with the help of scan- 
ning tunneling microscopy and scanning tunneling spec- 
troscopy, the presence of structure-dependent edge states 
of graphite can be observe d 16 ' 17 . As we know, some 
results of the edge states have been reported in the 
tight-binding model^Sr— , Dirac equation 2 ^— , the trans- 
fer matrix method (TMM) calculation^! - — , first princi- 
ples calculations^—, and most of them focus on the zero 
mode of the graphene nanoribbon, while some consider 
about the effects of disorder or the edge corrections. 

In this paper, we report a systematic study of the 
zero-energy and non-zero-energy edge states caused by 
the decoration on the edge by means of TMM. We an- 
alytically study the simple tight-binding model of semi- 
infinite ID lattice chain and the zigzag edged graphene 
(ZEG), and the consequent novel edge states that arise 
due to the edge decoration, which is the hopping inter- 
action on the edge being different from that in the bulk, 
caused as edge binding softening or stiffening. The ex- 
tended states and the localized edge states are also dis- 
cussed, and the analytical relations of the edge states 
with the edge decoration can be obtained as well. It 
might be applicable in future devices since novel edge 



states, particularly in the energy forbidden region, could 
significantly change the optical property, even the trans- 
portation. 

The rest of this paper is organized as follows. In Sec. 
IIIl we take a semi-infinite ID Peierls chain model as a 
fundamental case to show the relation of edge states with 
phase transition and describe how the analytical relations 
can be obtained efficiently through TMM for not only the 
extended states, but also the edge states. The sufficient 
and necessary conditions (SNC) for the existence of the 
edge states are determined. With the general discussion 
of the ID chain model, in Sec. IIII1 we map it to the ZEG, 
with applying the rigorous results in Sec. |TT]to obtain the 
edge states of ZEG with and without edge decoration. 
Finally, we give a brief conclusion in Sec. IIVI 



II. EDGE STATES AND PEIERLS PHASE 
TRANSITION 

As a typical example for studying edge states occurring 
in semi-infinite systems due to phase transition, we treat 
a semi-infinite ID lattice chain that can occur Peierls 
phase transition. It has been well-known that ID peri- 
odic atomic chain is not stable due to the electron-phonon 
interaction that drives chain to have pairing transition. 
In the case of half-filling electrons, the atoms in a periodic 
lattice chain like to pair together, forming Peierls pairing 
state. Its geometrical structure is shown in Fig. [TJ where 
one unit cell includes two atoms. In the tight-binding 
approximation, consider just the nearest-neighbor inter- 
action (n.n.), two different n.n. hopping constants are 
ti(= t ± A) and t2(= i T A) separately, where t is the 
hopping constant before Peierls phase transition, A(^ 0) 
describes the pairing phase transition just like the phase 
transition in ID molecule polyacetylene. A = when 
Peierls phase transition is suspended. In general, t\ ^ t% 
if A 7^ 0. The chain with Peierls pairing is equivalent to 
a ID model of atomic chain in which each cell has two 
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atoms, which are denoted in a cell as A and B. The 
edge of the semi-infinite chain is set at the left terminal. 
Without loss of generality, we assume that the lattice at 
left terminal is A denoted as 1, and label the unit cell 
with an integer number ordered from the left terminal 
to the right infinity. Since a non-Bravais lattice with two 
atoms per unit cell is bipartite, we have defined two kinds 
of fermion operators, ifiA and ifiB, where i is the inte- 
ger index labeling the unit cell. The Hamiltonian can be 
approximated by 

i^l i^2 

w 

where t ,ti,t 2 are three hopping parameters, and set all 
of them to be positive. The parameter to is introduced 
to describe edge decoration. From Eq. ([1]) , we can derive 
the dynamical equations 38 through defining an elemental 
excitation with energy E as follows 

{Eip 1A = t lfi 1B 

ElfiB = t (fiiA + h(f2A 

EifiA = hPiB + hfi-iB 

ElfitB = t 2 (fiiA + tilfii+iA (i > 2) 

According to above equations, we try to get the bulk 
wave-functions, and find out the physically meaningful 
localized states which decay from the surface. When 
E = or E 7^ 0, the above equations will be in differ- 
ent forms, so in the following subsection, we will discuss 
them separately. 
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(piA is the initial value of wave excitation with zero en- 
ergy at edge site. t 2 = t\ corresponds to the extended 
state in bulk. If t\ ^ t 2 , in order to obtain the physi- 
cally meaningful edge states, we must restrict t 2 jt\ < 1, 
so that when n goes to infinity, the wave-functions of 
sublattice A will approach to 0, showing the state is lo- 
calized near the edge region. Away from the surface, the 
amplitude of the wave-function will exponentially decay 
as e _n2a / AL . The localized length Al of the edge state 
is 2a/(lni 2 — lnti), where 2a is the periodic constant 
of the cell after Peierls pairing. Thus, we can conclude 
that zero-energy edge state can exist due to Peierls phase 
transition when t 2 < t\ (A ^ 0), and this result can 
be obtained by surface Green's function 3 ^. It is shown 
that no matter whether the edge decoration exists or not, 
zero-energy edge state is not affected by edge decora- 
tion £0(7^ t 2 ). ID Peierls pairing mode is an example to 
demonstrate that the studies of zero-energy edge state 
could show some evidence of the phase transition, even 
the mechanism of phase transition. 

In the following, we will study whether the non-zero- 
energy edge states can exist or not. According to Eq. 
@, after eliminate sublattice A and B separately, we 
can obtain 

E 2 lfiiA = t\ip 1A + tv,t\H>2A 

E 2 ip 2 A = {t\ + 1%)<P2A + tohlfiA + ht 2 (f ZA 
E 2 ipiA = [t\ + tQipiA + t\t^i-\A + tlt 2 ipi+lA,i > 3, 

(5) 
and 

{E 2 (fi 1B = (to + t\)(p 1B + tlt 2 lp2B, 

yE 2 Lp lB = (tf + t%)ipi B +tit 2 Lpi-iB +t 1 t 2 Lp i+1B ,i > 2 

(6) 

According to the first equation in Eq. we introduce 
a fictitious ipoA — into our discussion without lose of 
generality. Then we can turn Eq. ([5]) into matrix forms: 



FIG. 1: Schematic illustration of the lattice structure of ID 
semi-infinite chain with the closed (black) circles and the open 
(white) circles, consisting of sublattices A and B. 



Under the condition E = 0, 
the decoupled ones: 



Eq. ([2]) can be reduced to 



to<y5lB = 

toPlA +h(p2A = 

t 2 fiB+hfi-lB = 

h'PiA + tlfi+lA = 



(t>2) 
(i > 2) 



(3) 



From Eq. ([3]), it is easy to know that the wave- functions 
of all the type B sublattices will be equal to zero: {(fiiB — 
: i — 1,2, ••• ,00}; meanwhile the wave-function of 
sublattice A in the n unit cell will be: 



VnA = I - — 



to , 

"7" <PlA, 



n = 2,3,4' 



(4) 
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where T, T± and T 2 are defined as 

T = 



-1 




Ti = 



_*o 

t2 



-1 





(7) 



(8) 



with a = (E 2 



f2 



t 2 )/ht 2 



(E 2 - t 2 )/t h 



and j5 

Meanwhile, the wave-functions of the sublattices B sat- 
isfy the relations: 
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fi+lB 
ViB 



i-l I f2B 
<PlB 



V>2B 



E 
tt 



■flA 



(9) 



a > 2, we get < A_ < 1 and A + > 1, while A_ < -1 
and — 1 < A + < for the case of a < —2. We will discuss 
these two cases respectively. 

In the case of a > 2, it yields < A_ < 1, A + > 1. The 
necessary conditions to have a stable edge states localized 
at edge boundary must be T 2 — and 62 = that are 



and7 = (E 2 — t^— t\)/t\t 2 . The energy E must be real for 
the physical quantity, and all a, (3 and 7 are real numbers. 
Diagonalize the transfer matrix T to D = U~ 1 TU , and 
we have 



Vi+IA 
fiA 



UD i - 2 U- l T 1 T 2 




Vi+XB = UD i-i u -x 7 \e 



(10) 



with 



ViB 



D = 



Via- 



A_ 
A+ 



u- 1 = 



1 



Va 2 -4 



U = 



-1 
1 



A_ A+ 
1 1 



-A_ 



(11) 



(12) 



where A_ = |(a — \J o? — 4) and A + = |(a + V 'a 2 — 4). 
Eq. dTU|) can be rewritten as 



<fin+2B 



= u 



= u 



r x A!i N 




r 2A+ ^ 


VIA, 


9iA™ " 




e 2 A" 





(13) 
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Where I\ = C^ 1 (a/3 - |) + t/^ 1 /?, T 2 = 

f c/^ 1 /?, 6i - i/iiS + c/12 1 , e 2 = 

C/ 2 i 1 7 + t^22 1 - We must request 7^ 0, otherwise the 
wave function vanishes {fiA = ViB = : i — 1, 2, • • •, 00}. 

When |a| ^ 2, |A_| = |A + | = 1. It corresponds to 
the extended states in the bulk. Thus, all energy E 2 
satisfying the energy inequality \E 2 — t\ — t\\/\t\t 2 \ ^ 2 
must associate with the extended states. With the Eq. 
([5]) and Eq. ©, considering a continuous parameter k x 
to keep I cos(fc x a)| < 1, we can get the dispersion relation 
of the bulk states as: 



E 2 = tf + t\ + 2t 1 t 2 cos(k x a), k x a £ [-tt, tt]. (14) 

In the case, we have A± = e ±lfcxa . From Eq. (fT3|) . the k x 
belonging to K is clearly served as a wave vector of this 
extended mode propagating along the chain direction. 
The left edge makes the semi-infinite chain to have two 
waves interfered. One is the e jfcx ™ a , the other is e - lk * na . 

In the following, we focus on the solution of edge state 
when the condition \a\ > 2 is satisfied. In the case of 



[/-i L/3 _ *£j + u-10 = 

u£i + v£ = 

From Eq. (fTT]) . Eq. (Ti"2")) and the definitions of A±, above 
two equations can be changed to the following: 



(a + \f ct 



2„ 4 )/3-2^=0, 
4 - a + 2-i = 



(15) 



where a, (3 and 7 are the functions of energy E{^= 0). Eq. 
(fT5"j) is the necessary condition for the occurrence of non- 
zero-energy edge state. If Eq. (Ti~5l) can not be satisfied 
for any finite energy E, the non-zero-energy edge states 
could not exist. Inputting the definitions of a, f3 and 7 
into Eq. flT5"j) . we obtain 

(*o _ ^)(^ 2 - *o) = *o*i 

When the edge has no decoration, to = t 2 (y^ 0), there 
should be no finite solution of E for above equation since 
t\t\ 7^ 0. Thus we conclude that non-zero-energy edge 
state does not exist in ID Peierls chain without edge 
decoration. When to ^ t 2 , we have solution 



E 2 - t +t 2 t\/ [to - t 2 ) 



(16) 



Due to a > 2, we obtain E 2 > (t x + t 2 ) 2 > 0, which 
corresponds to a stable edge state. With the definition 
of 7 and a, we get 7 = a — (i§ — ^i) Ai*2- Taking it back 
to Eq. (irSI) . we have 



(t§ - t§) / tl t 2 = - 



> 1 



Thus the condition to get an additional stable non-zero- 
energy edge state is the inequality > t\ + t\t 2 . An 
appropriate manipulation on the edge decoration to can 
induce the non-zero-energy edge state. There is no con- 
straint showing the quantity relation between t\ and t 2 . 
From previous subsection, it has been known that the 
zero-energy edge state can exist with t 2 /ti < 1, so the 
zero mode and the non-zero-energy edge states can ap- 
pear at the same time. However, there will be non- 
zero-energy edge states only if t 2 ^ t\ together with 

ht 2 . 

Similarly we can get constraints for a < —2 



fj, ^ J-Z 

L ^ L 2 



(a - \Ju 2 -4)13 - 2 



\J a 2 



to 
h 

a — 27 
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They also yield the same energy condition: E 2 = ig + 
*o*i/ (*o — And besides 7 = a— (tp — t|) /tii2> when 
a < —2, it is easy to get: 



We can find that the non-zero-energy edge state can also 
exist under the condition: t\ < t 2 (t 2 — ti) that requests 
t 2 > t\ and results zero-energy edge state vanishing. 

From above two cases, we can see that non-zero-energy 
edge states are tightly related with edge decoration to and 
can be manipulated by decorating the edge atoms. 

We conclude that a zero-energy edge state will appear 
if Peierls transition happens (t^ < ti) regardless edge 
decoration at all. A non-zero-energy edge state can occur 
due to some edge decoration as following condition: 

E 2 = t 2 +t 2 t 2 /(t 2 -t 2 )(>0,t ^t 2 ) : 
t\ > t\ + tit 2 > 0; or tl < t 2 (t 2 - h), t 2 > h. 

(17) 

Non-zero-energy edge state can be manipulated with dec- 
orating the edge. 



III. EDGE STATES OF SEMI-INFINITE ZEG 

In this section, we would present some analytical re- 
sults for the edge states of ZEG. Although some results 
has been reported^ for uniform ZEG ribbon, and can 
be extended to semi-infinite ZEG through size scaling to 
the infinite width of ribbon, the effect of edge decora- 
tion on the edge states has not been reported. Thus we 
will supply some analytical results on the occurrence of 
non-zero-energy edge states and its energy dispersion in 
semi-infinite ZEG with edge decoration. In the previ- 
ous section, we have obtained some criteria on the ex- 
istence of edge states in semi-infinite ID Peierls model 
with and without edge decoration. In fact, some 2D and 
3D lattice models could be reduced into above ID model 
such as the model of ZEG to be discussed in this sec- 
tion. The geometrical structure of graphene is shown in 
Fig. [2 Without loss of generality, we denote the the 
first left line 1a as an edge, and take the lattice system 
in y direction to be infinite and periodic with constant 
a. Then the semi-infinite ZEG can be described by a 
line set with infinite number of periodic ID perpendicu- 
lar chain along y direction: {1a, 1b, n-A, ns, • • • }. The 
positions of sublattices A and B are labeled by two in- 
dices {i,Vi t A(B)} where i labels the number of perpen- 
dicular line from 1 to 00. After Fourier transformation 
for {yi.A(B)} with periodic constant a, in second quanti- 
zation representation, the model Hamiltonian with n.n. 
hopping can be expressed by a set of the Fermion oper- 
ators {(pi,A(B){k v ) ■ i = 1,2, ...,00}: 




FIG. 2: Schematic illustration of the lattice structure of semi- 
infinite zigzag edged graphene, the fulfilled circles represent 
the sublattice B, others are the sublattice A. The left first 
line 1a is the zigzag edge of the graphene. 



H = H t 4,A( k v) ( t>i-i,B{ky)+t4>l A (ky)(j)i,B(ky) 

k v ,i^2 

+ ^2*0^1 A {k y )(t>i,B{ky) + h.c, (18) 

Ky 

where k y is a good quantum number, t = 2t cos(k y a/2) 
and to — 2t' cos(k y a/2). We have introduced a simplified 
edge hopping decoration t' at first line. When t' = t, it 
becomes an ideal semi-infinite ZEG without edge deco- 
ration. In general, t' 7^ t that describes some edge deco- 
ration, and we set all the hopping interaction positive as 
well. The wave function at site of A(B) sublattice 

is written as ^(B)( k y) = ex P( i k y yj,A(B))<l>i,A(B)(k y )- 

For each fixed k y , the Hamiltonian of Eq. (ITS)) is equiv- 
alent to the ID Peierls chain in Sec. [TTI when we denote 
to = to, ti = t and t 2 — t. Thus the conclusion in Sec. |TT] 
can be applied to study the ZEG. 

Comparing with standard energy spectrum of bulk 

1 E 2 -4t 2 cos 2 k v a/2-t 2 , 2 ^ a 

graphene, a = OT , L , when a ^ 4, we can 

G " ' 2t z cos fc y a/2 5 ^ 1 

get the dispersion relation for extended states: 

E 2 = tl+t 2 2 + 2t 1 t 2 cos{Vik x a/2) 

= t 2 [l + 4cos 2 (k y a/2) + Acos(k y a/2) cos(V3k x a/2)} 

(19) 

This is the standard formula of the energy spectrum of 
the infinite bulk graphene. 

In the following, we would focus on the edge states. 
From Sec. [TTI the necessary condition for the existence 
of zero-energy edge state is | *2 / * 1 1 < 1 that directly re- 
sults the condition |2 cos(fcj,a/2)| < 1 for ZEG. It clearly 
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the non-zero-energy edge states depends on the decora- 
tion on the edge, the ratio between the edge hopping and 
the bulk hopping should be within a certain threshold, so 
that we can detect the finite-energy edge states. Compar- 
ing the results obtained above with those by the means of 
the exact diagonal method on a finite-size zigzag nanorib- 
bon model, the results are accordance with each other, 
as shown in Fig. [3] 



FIG. 3: (color online) By exact diagonalization method, the 
dispersion relation of the bulk states and the edge state are 
depicted (blue line), while the dispersion relation of the edge 
state obtained through TMM is depicted in the dashed line 
(red dash-line). Parameters: (a) t' — 0.5t, (b) t' — 2.0t. 



shows the well-known conclusion that a zero-energy edge 
state exists in the region k y a <E (27r/3, 4tt/3) for ZEG. 
But the point of k y a = ±tt must be excluded since 
t% = 2tcos(k y a/2) — 2icos(7r/2) = that corresponds 
to an internal mode within the isolated chains. From the 
conclusion of ID model in SecJITl we also can conclude 
that the zero-energy edge state is stable. 

Just as the discussion shown in the above section, the 
SNC for the existence of the non-zero-energy edge states 
should be discussed with a 2 > 4. With similar discus- 
sions, we can get the dispersion of energy spectrum of 
non-zero-energy edge states as follows: 



E 2 = At'n 



(k y a/2) 



t' H 



(20) 



while the SNC for the existence of those states must sat- 
isfy: 

t' < t, < 2^2%^ < I cos (k y a/2)\ «C 1 

=>t 2 >2t' 2 ; 
or (21) 
t' > t, < 2(t /_ t 2) < I cos {k y a/2)\ < 1 
=>• t 2 < 2t' 2 /3 

As above shown, we can easily see that the existence of 



IV. CONCLUSIONS 

We have studied the semi-infinite ID lattice chain with 
Peierls phase transition and ZEG model respectively, 
with the help of TMM. We focus on the zero-energy 
edge states and non-zero-energy edge states, under the 
condition of the edge decoration, which is the electron 
hopping energy being different from that in bulk. In the 
ID model, we conclude that a zero-energy edge state will 
appear if Peierls transition happens regardless edge dec- 
oration at all. A non-zero-energy edge state can occur 
due to the edge decoration. Then we map the ID case 
into the ZEG. Suppose the electron hopping energy is 
£q, different from that t in bulk because of the boundary 
effects. It is found that the existence of the zero-energy 
edge states is independent of the magnitude of t' . In ad- 
dition to the zero-energy edge states, there can exist edge 
states of finite energy if t > y/2t' or t < y/2/3t' . These 
edge states appear as pairs with energy ±E outside of 
the continuum conduction and valence bands. 
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